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In the framework of inflationary models with non-vacuum initial states for cosmo- 
logical perturbations, we study non-Gaussian signatures on the cosmic microwave back- 
ground (CMB) radiation produced by a broken-scale-invariant model which incorporates 
a feature at a privileged scale in the primordial power spectrum. 
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1. Introduction 

The common belief that the CMB is Gaussian distributed can be directly traced 
back to the generic assumption that the quantum fluctuations of the inflaton field 
are placed in the vacuum state Relaxing this assumption might lead to detectable 
signatures in various astrophysical tests, most interestingly in future CMB and 
large-scale structure observations. In this note, we study CMB non-Gaussian signa- 
tures predicted within inflationary models with non-vacuum initial states for cos- 
mological perturbations. The model incorporates a privileged scale, which implies 
the existence of a feature in the primordial power spectrum. The model predicts a 
vanishing three-point correlation function for the CMB temperature anisotropics^. 
We here focus on the first non- vanishing moment, the CMB four-point function at 
zero lag, namely the kurtosis, and compute its expected value for different locations 
of the primordial feature in the spectru m, a s suggested in the literature to conform 
to observations of large scale structure 



2. Two-point correlation function for non-vacuum initial states 

We consider non- vacuum states for the cosmological perturbations of quantum me- 
chanical origin. Let 2?(<r) be a domain in momentum space, such that if k is between 
and a, the domain T>(<r) is filled by n quanta, while otherwise T> contains nothing. 
The state |\&i(<7, n)) is defined by 

(J )n 

\*i(<r,n))= J] M^flOO ® I°p>= ® ' nk > ® I°p>- f 1 ) 

ker>( CT ) pg£>(<r) keX>(er) P^V{ct) 
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The state |nk) is an n-particle state satisfying, at conformal time r\ — c^n^) = 
y^Kn — l)k) and c k |rik) = \/n + l|(n + l)k)- We have the following property 

{*i(<T,n)\*i((j',n'))=5(a-(j')5 nn >. (2) 

It is clear from the definition of the state l^i) that the transition between the 
empty and the filled modes is sharp. In order to "smooth out" the state we 
consider a state |*2) as a quantum superposition of In doing so, we introduce 
an, a priori, arbitrary function g(a; fcb) of a. The definition of the state \ 1 t l 2( n , &b)) 

r+oo 

|tf 2 (n,fcb))=/ do- fl (a; fcb)|*i (o-.n)), (3) 
Jo 

where g{a;k\ > ) is a given function which defines the privileged scale fcb- We as- 
sume that the state is normalized and therefore J Q °° g 2 (a; kb)da = 1. In the state 
\^i(a,n)), for any domain V one has 

(tfiforOlcpCqltfifon)) = n)|ct C t |*i (tr,n)) = 0, (4) 

(*i(tr,n)|cpct|*i(a,n)) =n5(qeX>)5(p-q)+J(p-q), (5) 
(*i(tr,n)|ct Cq |*i(<r,n))=n5(qeI>)5(p-q). (6) 

In these formulas, <5(q € I?) is a function that is equal to 1 if q £ T> and otherwise. 
These relations will be employed in the sequel for the computation of the CMB 
temperature anisotropies for the different non- vacuum initial states. 



2.1. Two-point function of the CMB temperature anisotropy 

The spherical harmonic expansion of the cosmic microwave background temperature 
anisotropy, as a function of angular position, is given by 

St f St 

-jrie) =J2a em Y em (e) with a lm = J dQ. e — (e)Yf m (e). (7) 

As we are interested in a non-Gaussian signature of primordial origin we will be 
focusing on large angular scales, for which the main contribution to the temperature 
anisotropy is given by the Sachs- Wolfe effect, namely, 5T/T(e) ~ (l/3)$[r/i ss , e(?7o — 
f]\ss)], where x) is the Bardeen potential, while rj and ?7i ss denote respectively 
the conformal times now and at the last scattering surface. Note that the previous 
expression is only valid for the standard Cold Dark Matter model (sCDM). In 
general, we might also be interested in the case where a cosmological constant is 
present (ACDM model) since this seems to be favored by recent observations. Then, 
the integrated Sachs- Wolfe effect plays a non-negligible role on large scales and the 
expression giving the temperature fluctuations is not as simple as the previous one. 

In the theory of cosmological perturbations of quantum mechanical origin, the 
Bardeen variable becomes an operator, and its expression can be written as 



c k (r?i)/ fc (r?)e ikx + 4(r ?i )/ fc *( ?? ) e - tkx 



(8) 
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where £p\ = (Gh) 1 / 2 is the Planck length. In the following, we will consider the 
class of models of power-law inflation since the power spectrum of the fluctuations 
is then explicitly known. In this case, the scale factor reads a{rf) = lo\f]\ X ^^ ■> where 
(3 < —2 is a priori a free parameter. However, in order to obtain an almost scale- 
invariant spectrum, f3 should be close to —2. In the previous expression of the scale 
factor, the quantity £q has the dimension of a length and is equal to the Hubble 
radius during inflation if j3 = — 2. The parameter l§ also appears in Eq. JHJl. The 
mode function fk(r]) of the Bardeen operator is related to the mode function fJ>k(v) 
of the perturbed inflaton through the perturbed Einstein equations. In the case of 
power-law inflation and in the long wavelength limit, the function fk{vi) is given 
in terms of the amplitude A s and the spectral index n s of the induced density 
perturbations by 



fc 3 |/ fc | 2 = A s k r ' 



(9) 



Using the Rayleigh equation and the completeness relation for the spherical har- 
monics and after some algebra we get 



atn 



^e M ' 2 I dk 



At this point we need to somehow restrict the shape of the domain T>. We assume 
that the domain only restricts the modulus of the vectors, while it does not act on 
their direction. Then, from Eq. 1)10(1 . one deduces 



(*i(cr, n) 



eimia| 2m3 |*i(o",n)) = 



Ct x +2nD { e l\a) 



dk 

f £ [k( m -m ss )]k 3 \f k \ 2 T 



' _ " i a L m 2 i (11) 

(12) 



Thus, the multipole moments Cl , in the state are given by C\ L> [a 



(i)/ 



2nDg 1 \a) , where Cg is the "standard" angular power spectrum, i.e., the multipole 
obtained in the case where the quantum state is the vacuum, i.e., n = 0. Let us 
calculate the same quantity in the state l^)- Performing a similar analysis as the 
above one, we find 



f 2 



(*2(n,fcb)|a^ imi al m2 |* 2 (n,fcb)) = - w 



D 



(2) 



7T 



Ct 



2nD 



(2) 



-\-oc 



J 2 +1/2 (k)h(k)k ns ~ 3 dk 



e 1 e 2 dm 1 m 2 ! (13) 
(14) 



where, to reach this eqn, we defined g 2 (a;k\>) = Ah/ 'da [we will see below that 
this function h(ky,) cannot be arbitrary] and we integrated by parts, and then we 
defined h(k) = /i(oo)[l — h(k)/h(oa)]. In this we have not assumed anything on 
h(oo) or h(0). We see that the relation g 2 {k) = dh/dk requires the function h(k) to 
be monotonically increasing with k. It is interesting that, already at this stage of 
the calculations, very stringent conditions are required on the function h(k) which 
is therefore not arbitrary. This implies that the function h(k) which appears in the 
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correction to the multipole moments is always positive, vanishes at infinity and is 
monotonically decreasing with k. An explicit profile for h(k) is given in Fig.^ The 
total power spectrum of the Bardeen potential can be written as 

h(k) 



fc d |$ fc | « A s kns ~ 1 1 + 2n/i(oo 



1 



h(oo) 



A s k n s~ 1 [l + 2nh(k)}. (15) 



Observations indicate that n s ~ 1 and for simplicity we will take n s = 1. As we 

have seen previously, we can write the multipole moments in the state as 

(2) (2) 

C) = Ce + 2nD) . Substituting the well-known expression for the Ci's and the 



definition of D\ given by Eq. |Q , one finds that the coefficients C\ are given 

by 



C, 



(2) 



1 r(3-n s )r[£+(n s -l)/2] 
2 1 2 3 -" s T 2 [(4 - n s )/2]T[£ - (n s - 5)/2] 



2nD 



(2; 



(16) 



As a next step, one has to normalize the spectrum (need to determine the value of 
A s ). We choose to use the value of Q rms _ PS = T [5Cf ] /{^)] 1/2 {l P \/£ ) ~ !8fiK 
with T = 2.7K measured by the COBE satellite. Thus, we compute the quadrupole 
and then 

-1 -1 



8Q 



rms-PS c 



T 2 



i 2 



1 

(ill 



2nD. 



(2) 



for ri s = 1. The band power STg gives 



Qrms-PS 



12 



l + 2mr£{£+l)D 



(2) 



1 + 12nnD. 



(2) 



(17) 



(18) 



The n-dependence in the above expression is the correction due to the non- vacuum 
initial state. We easily check that if n = the corresponding band powers are 
constant at large angular scales. 

Finally, we calculate the two-point correlation function at zero lag in the state 
1^2)- Using Eqs. Q), i|13|) . the second moment, /i2, of the distribution is given by 

1 2\ 



M2 



ST. , 
jr(e) 



= 72" 2^ 



4tt 



(19) 



Once we have reached this point, an obvious first thing to do is to check that the two- 
point correlation function calculated above is consistent with present observations. 



2.2. Comparison with observations 

Among the available observations that one can use to check the predictions of 
theoretical models, two are key in cosmology: the CMB anisotropy and the matter- 
density power spectra. We will not study in details all the predictions that can 
be done from the two-point correlation function since our main purpose in this 
note is to calculate the non-Gaussianity which is a clear specific signature of a non 
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Fig. 1. The function h(k) for different values of a (left panel). On the right panel we show 
the matter power spectrum normalized to COBE for different numbers n of quanta in the initial 
state. The cosmological parameters are those corresponding to the sCDM model, namely, h = 0.65, 
Q\ = 0, Cb = 0.05, f2 c( j m = 0.95 and n s = 1. The parameters describing the non vacuum state 
are k^ hyB = 0.052/iMpc" 1 and a = 2.5. The data points represent the power spectrum measured 
by the PSCz survey. 



vacuum state. So, we just compute the matter power spectrum to demonstrate that 
it fits reasonably well the available astrophysical observations for some values of the 
free parameters. In addition, this illustrates well the fact that, using the available 
observations, we can already put some constraints on the free parameters. A simple 
ansatz for the function h{k) is represented in Fig. ^ and can be expressed as 



1 



, / , k 
tanh |am — 

k h 



(20) 



With it, the matter power spectrum today, after taking into account the transfer 
function T(k) which describes the evolution of the Fourier modes inside the horizon, 
can be written 



P{k) = T 2 {k) 



ms-PS 



T 2 



— + 2nD ( i\k h ) 

07T 



1 + 2nh(k) 



phys 



The sCDM transfer function is given approximatively by the numerical fit 



5 



T(k) 



ln(l + 2.34g) 
2Mq 



1-1/4 



' + (16. lq) 2 + (5.46g) 3 + (6.71g)' 



where q and the shape parameter T can be written as 

q = fc/KM^Mpc' 1 ] r = nohe-^-^^/^ 



(21) 



(22) 



(23) 



where fin. is the total energy density to critical energy density ratio and fib rep- 
resents the baryon contribution. Or, more explicitly, we take £Iq = fl\ + O m = 
SIa + ficdm + f^b- We have now normalized the matter power spectrum to COBE. 
It is important to realize that the above procedure only works for the sCDM model 
since we have used the Sachs- Wolfe equation. The sCDM matter power spectrum 
is depicted in Fig. ^ The measured power spectrum of the IRAS Point Source 
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Catalogue Redshift Survey (PSCz) I™ has also been displayed for comparison. One 
notices that the effect of the step in h(k) is to reduce the power at small scales 
which precisely improves the agreement between the theoretical curves for n ^ 
and the data. Let us remind at this point that the shape of the function h(k) has 
not been designed for this purpose and comes from different (theoretical) reasons. 
Therefore, it is quite interesting to see that the power spectrum obtained from our 
ansatz fits reasonably well the data. This plot also confirms the result of Ref. 
namely that the number of quanta must be such that (for sCDM) 1 < n < 10, i.e., 
it cannot be too large. 



3. Four-point correlation function for non-vacuum initial states 



It is time to proceed with the calculation of the four-point correlation function. 
Briefly one has to first perform the calculation for the state l^i), and thengener- 
alize it for |\I/ 2 ). After some lengthy but straightforward algebra, one finds ^ 



|* 2 (n,fc b )) 



pi 



C tl C ta + 2nC tl Df] + 2nC e2 D<£> + 4n 2 F$ 2 



+ (— l) mi+ ™ 2 
_|_/_^\mi+m 3 

Z 'H' J ■+■ L )- Cj £ 1 e 2 e 3 ei ri e 1 £2e 3 i 4 e 



C ei C i2 + 2nQ 1 £)g ) + 2nC l2 D^ + 4n 2 F$ 2 
C tl C i3 + 2nC ei D<£> + 2nC e3 D^ + 4n 2 F$ 3 



$£\£i & £2^3 , — m&$m2 ,— 7713 



5£ 1 £ 2 5£ 3 £ i 6 mi} - m2 8 m3y - m4 



! + (_!)/!+<» + (_!) 



(24) 



? (2) 

tlt2 



+ °° _ d 

dah(a) — 
da 



E 



(2) 



kAKvo - %*)] ■ ■■jeAHvo - mss)]h(k)k 3 \f k \ 



k dk 
T 



(25) 
(26) 



We are now in a position to calculate the CMB excess kurtosis. In order to establish 
an analytical formula for it, one just needs to use the equation linking ag m and 
5T/T and play with the properties of the spherical harmonics. Explicitly, the excess 
kurtosis is defined as K, = /i4 — 3/i| , where the second moment has already been 
introduced and where the fourth moment, /i4, of the distribution is defined as 
/i4 = {K) with K = [^(e)] 4 . An important shortcoming of the previous definition 
is that the value of K, depends on the normalization. It is much more convenient 
to work with a normalized (dimcnsionless) quantity. Therefore, wc also define the 
normalized excess kurtosis as 



2^ = ^-3 



(27) 



which is the one more commonly used in the literature. We have performed the full 
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Fig. 2. The normalized excess kurtosis parameter Q in terms of the privileged (comoving) wave- 
number fct, for a particular representative set of parameters: a = 5 for the sharpness of the weight 
function signaling the privileged scale, and n = 2 quanta in the non-vacuum initial state for the 
cosmological perturbations. 



numerical resolution of the normalized excess kurtosis parameter which we show in 
Fig. [3 As we see from it, Q ~ —4.24 x 10~ 4 is an asymptotic value, provided we 
concentrate on the middle and big values of the built-in scale fcb- An analytical order 
of magnitude estimate of this can be seen in Ref. "4 The fact that the numerical 
estimate does not depend on k\, is confirmed by the plot, except for small values 
of the wave- numbers. In fact, this shows that the quantity Q does not depend very 
much on the free parameters. One can see this is true for parameters a and n since, 
using the above mentioned analytical estimate, we find, for n = 1, Q ~ — 3.77x 10 -4 
and for n — > oo, Q ~ —4.25 x 10 -4 . Since we know that this result does not depend 
on the details of the weight function h(k), we conclude that the asymptotic value 
obtained above is a generic value, at least for large values of fch>. I n particular, this 
is true for kb ~ 300 which corresponds to the built-in scale located roughly at the 
privileged scale in the matter power spectrum selected by the redshift surveys of 
Ref. El. Another important remark is that the excess kurtosis is found to be negative. 



4. Discussion 

The above paragraphs show that the introduction of non- vacuum initial states leads 
to a clear prediction for the excess kurtosis parameter. However, this by itself does 
not imply that the signature will be observable. The signal to theoretical noise ratio 
has to be considered or, stating it differently, we need to make the comparison of 
the signal against the cosmic variance. The cosmic variance quantifies the theo- 
retical error coming from the fact that, in cosmology, observers have only access 
to one realization of the ST/T stochastic process whereas theoretical predictions 
are expressed through ensemble averages. The specific computation of the cosmic 
variance for the model under consideration was performed'" and found to be some 
four orders of magnitude higher than the signal itself. This is no surprise: we are 
interested in a non-Gaussian signature of primordial origin, and thus we are focus- 
ing on large angular scales, for which the theoretical uncertainties are the highest. 
Shifting to intermediate angular scales, a stronger signal would be obtained; how- 
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ever, in that case secondary sources would be more difficult to subtract and thus 
the transparency of the effect would be compromised. 

In sum, the excess kurtosis is found to be negative and the signal to noise ratio 
for the dimensionless excess kurtosis parameter is equal to \ S/N\ ~ 4 x 10 -4 , almost 
independently of the free parameters of the model. This signature turns out to be 
undetectable and therefore we conclude that, subject to current tests, Gaussianity 
is a generic property of single field inflationary models. 
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